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<^ ; Abstract 

^ . We revisit a model of the two-scalar system proposed previously for un- 

derstanding a small but nonzero cosmological constant. The model pro- 
vides solutions of the scalar-fields energy ps which behaves truly constant 
^ ' for a limited time interval rather than in the way of tracker- or scaling- 

O ' type variations. This causes a mini-inflation, as indicated by recent ob- 

^ . servations. As another novel feature, ps and the ordinary matter density 

On ! pm fall off always side by side, but interlacing, also like (time)~^ as an 

overall behavior in conformity with the scenario of a decaying cosmolog- 
ical constant. A mini- inflation occurs whenever ps overtakes pm, which 
may happen more than once, shedding a new light on the coincidence 
W). problem. We present a new example of the solution, and offer an intu- 

itive interpretation of the mechanism of the nonlinear dynamics. We also 
discuss a chaos-like nature of the solution. 



1 Introduction 

Sometime ago we proposed a theoretical model for the system of two scalar fields as 
an effective way to understand a small but nonzero cosmological constant, as indicated 
by a number of observations. The indication is now even stronger in particular due to 
the recent results on type la supernovae at high redshift In accordance with this 

development we provide in this paper, (i) an example of the cosmological solution suited 
for the recent data, (ii) an intuitive interpretation of the underlying mechanism of the 
model, (iii) a novel way of viewing the "cosmological coincidence problem." We also add 
comments on the related issues including a generalization, chaos- like nature of the solution 
and choice of a physical conformal frame. 

It seems useful to outline beforehand how we arrived at this particular model, which 
was conceived in two steps. We first started with the prototype Brans-Dicke (BD) model 
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with a constant A added ||^. We apply a conformal transformation (Weyl rescaling) to 
remove the non-minimal coupling term. We say we have moved from a conformal frame 
(CF) called J frame after Jordan to another called E frame after Einstein (denoted by * 
hereafter if necessary). In the latter CF we have an exponential potential 

V{a) = Ae-^'^'^, (1) 
where a is a canonical scalar field related to the original field by 

= r'/'e^^ (2) 
The constant ^ is related to the original definition a;by4^c<j = e = ±l, and 

r' = 6 + er\ (3) 

which should be chosen to be positive so that a is a (non-ghost) normal field. 

This exponential potential serves to implement the scenario of a "decaying cosmo- 
logical constant," Agff = Po- ~ (with the cosmic time in E frame), allowing us to 
understand why today's A is smaller than the theoretically natural value by as much as 
120 orders of magnitude, where p„ is the energy density of cr 0. 

We find, however, that pa which falls off in the same way as the ordinary matter density 
p^m is not quahfied to explain what the recent observations indicate. We need p„ falling 
off more slowly than or more preferably staying constant, at least for some duration 
of time. A scalar field expected to show this feature, now often called quintessence, has 
been a focus of many studies ||7|-|T0|. 

In this connection we notice that the (analytic) solution resulting in p^ ~ is 
an attractor realized asymptotically. We know, on the other hand, that there is an 
interesting transient solution, which allows a to stay nearly constant temporarily. This 
is a consequence of a fast decreasing potential and the frictional force supplied by the 
cosmological expansion 0,^^ . The plateau behavior of a lasts until comes down to 
be comparable with the overdamped Po-, when the system goes into the asymptotic state. 

It may appear that this plateau mimics a cosmological constant. It does, but not to the 
extent required to fit the observations, because a is not persevering enough to keep staying 
constant as it is seduced to move into the asymptotic phase too early. This is a place 
where we enter the second step, in which we expect occasional "small" deviations from 
the dominant behavior ~ a crucial pattern to maintain the decaying cosmological 
constant scenario. 
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After a desperate search for the success, we decided to introduce another scalar field 
$ that couples to a through the potential 

= e"^'^"y((T,<l>), with \/ = A + im2$2 [1 +7sin(Ka)] , (4) 

with m,7 and k constants [||, as also shown graphically in Fig. |l]. This is supposed to 
be a potential in E frame. For $ = 0, we go back to the exponential potential ([^) derived 
from A simply by a conformal transformation. At this moment, we have no known theory 
at a more fundamental level from which $ and V^(cr, $) derive, nor we claim that is a 
final result. We nevertheless emphasize that we then find solutions which are qualitatively 
different from any of the results known so far, expecting to open up a new perspective on 
the nature of the cosmological constant problem. 

In Section 1 we describe the model very briefly. We then show a typical example of the 
solutions in Section 2, comparing the results with observations, before we go into some 
detailed discussion of the underlying mechanism in Section 3. Final Section 4 is devoted 
to other related discussions. 




2 The model 

We reproduce the basic field equations in Refs. [||]. In terms of r = Int* and 6 = Ina* 
with a^, the scale factor in E frame, the cosmological equations in spatially fiat Robertson- 
Walker universe can be written as 

U>' + V{a,<^\ (5) 
= 0, (6) 

uo I 

$" + (36' - 1)$' + ^'e-^^'^^ = 0, (7) 
pL + 46'p,™ = 0, (8) 

where a prime implies a derivative with respect to r. We seek solutions of a spatially 

uniform a. We also assumed the radiation-dominated universe, for simplicity at this 

moment. Due to the relation H^, = {da^/ dt^) / a^, = t~^b', we find that choosing r as an 

independent time variable has an advantage that the coefficients of the frictional force in 

(^) and are now constant as far as a^,(t*) expands according to a power-law. On the 

other hand, the potential in these equations is multiplied by tl = e^'^. We use the reduced 

Planckian unit system.^] 

^Units of length, time and energy are 8.09 x lO^'^'^cm, 2.70 x 10^'^^s and 2.44 x lO^^GeV, respectively. 
Note that the present age of the universe ^ 10^°y corresponds to ^ 10^". 
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We choose the initial time of integration to be 10^'', somewhat later than the time 
when the reheating process is supposed be completed, though the true initial conditions 
for the classical fields should be given at a much earlier time. 

3 An example 

We give an example of the solution of (|^)-(|§), as shown in Fig. ^. We came across this 
solution by a somewhat random search for the parameters and the initial values. We 
have not attempted to scan the whole parameter space, though the values have been 
constrained essentially of the order one in the reduced Planckian unit system. 

At the same time, we fine-tuned the parameters and initial values moderately in order 
to obtain considerable amount of Q\ at the age around 10 Gy, where we find a crossing 
between ps and p*m. This naturally induces an extra acceleration of the scale factor a*, 
a "mini-inflation." In Fig. ^, a magnified view of the lower diagram of Fig. ^ we present 
more detailed behaviors of the densities, from which we can calculate and h, the 
Hubble constant in units of 100 km/sec/Mpc, for several values of the assumed present 
age of the universe, as listed in the "first set" of Table |1|. 

It may appear that Qa and h shown in the first set are shifted slightly to the higher 
side of those of the recent data; ~ 0.7, h = 0.65, t^Q = 1.3 x 10^°y as summarized in 
Instead of searching for another set of parameters, we exploit the fact that the eqs. (|§)-(|§) 
are invariant under the rescaling t^, — > Xt^, combined with A — > A^A, m —> Am, p^m.s 
X^P*m,s- By choosing A = 0.7762 to shift logt* = 60.11 to 60.00, for example, we obtain 
the "second set" of Table yielding somewhat lower values of VL\ and h. It seems likely 
that we are able to adjust ourselves to almost any of the observational results which will 
be updated in the future with less uncertainties. 



t*o(10^V) logt*o 


First set 
i^A h 


Second set 

r^A h 


1.1 60.11 

1.2 60.15 

1.3 60.18 

1.4 60.22 


0.62 0.81 
0.67 0.77 
0.72 0.74 
0.76 0.72 


0.46 0.73 
0.52 0.69 
0.56 0.65 
0.62 0.63 



Table 1: The assumed ages t^o, logt*o in the Planckian unit system, Qa and h for the 
solution in Figs. ^ and ^ (First set), and for the rescaled age and other parameters 
(Second set). 
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We list other marked features shown in Fig. 0. (i) We have two plateau periods and 
the associated two mini-inflations. We have other examples with even more mini-inflations 
(with none as well). This would imply that what we are witnessing at the present time 
may not be a once-for-all event, but can be only one of the repeated phenomena in the 
whole history of the universe. In this way we can make the coincidence problem less 
severe, (ii) The mini-inflation which has just started at the present time will not last 
forever, (iii) In lower diagram of Fig. 0, we find that around the time of nucleosynthesis 
(logt* ~ 45), ps stays several orders smaller than p^^n- Corresponding to this, the middle 
diagram of Fig. ^ shows that the effective exponent a* = t^,H^, = b' stays at 0.5. It thus 
follows that the conventional analysis of nucleosynthesis is left unaffected by the presence 
of the scalar fields. This is a crucial criterion by which we can select acceptable solutions, 
(iv) Around the present epoch, our solution does not behave like a scaling or a tracker 
solution [0. Our solutions correspond exactly to the analyses in Refs. 0,0] including 
simply a cosmological constant. The proposed tests on the equation of state p/ p might 
be interpreted as selecting models of the scalar field(s) rather than distinguishing them 
from the models of a purely constant A. (v) A sufficiently long plateau of a may explain 
why the time-variability of some coupling constants as suspected from the time evolution 



of a have been constrained to the level much below t^^ ~ 10 ^'^y ^ [p!^,p 



4 The mechanism 

Let us recall that already in the system of a alone there is a tendency that a is likely 
overdamped, because once it begins to move toward infinity, the exponential potential 
decreases so fast that it is only decelerated by the cosmological frictional force eventually 
to a standstill on the middle of the potential slope. With $ 7^ 0, on the other hand, there 
is now a chance that a is trapped in one of the minima of the sinusoidal potential. With 
a still persistent frictional force, $ will stay also perching on the slope. This will continue 
to provide the potential wall of a, thus enhancing the chance for a to stay sufficiently 
long, giving a sufficiently large fl\. 

In the mean time, the strength of the potential in (|^) will grow due to the factor t^. 
The potential wall gets increasingly steeper, helping the energy of a to build up. For the 
same reason the force acting on <I> grows also. The exponent /(r, a) = 2T — A(a, which has 
been negative for a which had advanced sufficiently, will increase again. If f{T,a) ~ is 
reached, $ is pushed downward for the central valley $ = 0. Consequently, the potential 
wall which has been confining a finally collapses, releasing the accumulated energy to 
"catapult" a forward, as shown in Fig. ^ a magnified plot of the upper diagram of Fig. 
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B 

We may expect that a is again decelerated and trapped by a potential minimum, which 
has in fact quite wide a basin of attraction, hence repeating nearly the same pattern as 
before, as shown in Fig. |^. It may also happen, on the other hand, that a fails to be 
trapped. Then without a potential wall for forced confinement, a will goes smoothly into 
the asymptotic stage as in the model with a single scalar field. Figs. provide an 
example of the behavior of this type. We find no crossing between the two densities, thus 
with VL^ reaching only to 0.29, hence without extra acceleration of the scale factor around 
the present epoch. 

Since the probability of being trapped is limited, however, the system will be, perhaps 
after some number of repeated mini-inflations, ultimately in the asymptotic behavior for 
any initial values. 

5 Discussions 

As we admitted, the potential (^) discovered on a try-and-error basis is still tentative. A 
crucial ingredient is obviously the presence of a sufficient number of minima with respect 
to 0", each with a considerable basin of attraction. These traps should also depend on 
another dynamical field. Other more general potentials might be found resulting in the 
unique feature that mini-inflations can occur repeatedly, in connection with the interlacing 
pattern of ps and p^rn- 

There are several other points to be discussed. First, there is an issue on what CF 
is a physical CF. Our calculations have been carried out in E frame. In this connection 
we point out that in the prototype BD model the scalar field is assumed to be decoupled 
from the matter part of the Lagrangian in J frame in order to maintain Weak Equivalence 
Principle (WEP). It then follows that masses of matter particles are time-independent 
even if the scalar field evolves. In this sense J frame should be a physical CF, because 
we analyze nucleosynthesis, for example, on the basis of quantum mechanics with particle 
masses taken as constant. We find, however, that the prototype BD model with A included 
entails an asymptotically static (radiation-dominated) universe, to which the universe 
approaches oscillating or contracting 0]. Although this conclusion depends on the simplest 
choice of the non-minimal coupling term as well as a purely constant A, we must modify 
the model in a rather contrived manner .0 

As another approach we proposed to modify the mass terms of matter fields guided 
by the invariance under global scale transformation (dilatation) [Q, making E frame an 

^If wc multiply A by 0^ {i ^2), the scale factor a{t) in J frame is found to behave like a{t) = with 
a={l/2){l/{t-2)). 
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approximately physical CF. WEP is then violated but can be tolerated within the con- 
straints from the fifth- force- type experiments. Leaving the details in Ref. ||^, however, 
we emphasize that we can evade large part of the complications in the choice of CF if we 
choose the type of solutions as exemplified in Fig. 0, because a stays constant over the 
period which covers nucleosynthesis and the present time. Not only o and hence particle 
masses at the time of nucleosynthesis remain the same as today, but also a constant scalar 
field makes any difference among different CFs trivial. 

Strictly speaking, however, we must include dust matter to solve the equations beyond 
the time of recombination, t^^eq ~ 10^^. For this purpose we in fact added Cd^*P*m and 
— Cd^*o"'p*m to the right-hand sides of (|]) and (||), respectively, simultaneously replacing 46' 
on the left-hand side of (||) by 36'. We used the coupling strength estimated tentatively as 
Cd = 0.005 determined as ~ 1/4 of the coefficient ~ 0.02 in eq. (58) of Ref. 0. Effects 
of these modifications turned out rather minor, though we recognize a slight distortion of 
the otherwise fiat ps for 50<, logt* <,57 in the lower diagram of Fig. ^ We consider this 
as the maximally expected effect. 

As another aspect of our solution, we should mention the sensitivity on initial values. 
The absence of the second mini-inflation in Fig. ^ is a consequence of a slight change 
of the initial value 0\ from 6.7544 to 6.7744 chosen in Fig. |^. We have in fact the 
solutions without any mini-inflation for in between. According to our experience, we 
come across a solution of two mini-inflations rather easily if we scan a\ over the range 
Ao"i ~ 0.6 (corresponding to the change 27r of Ka\ with a pitch 0.01, for example. More 
details will be reported elsewhere. 

Traditionally we have made it a goal to make predictions depending as little as possible 
on as few as possible number of initial conditions. Our solutions do not seem to meet this 
condition. We have at least four initial values for a and Moreover, the result depends 
on them very sensitively. One may suspect chaotic behaviors in nonlinear dynamics of 
more than three degrees of freedom; we have at least four of them from the two scalar 
flelds. 

As it turns out, however, the flnal destination of the solutions in phase space is a 
flxed-point attractor instead of any kind of strange attractors, as long as the cosmological 
frictional force continues to be effective. Our solutions show only chaos- /2A;e behaviors 
rather than authentic chaotic behaviors as discussed in Refs. [1^. A conventional de- 
terministic view may not be maintained. In contrast to the long-held belief in physics, 
it seems that we are forced to accept a new natural attitude if the nonlinearity of the 
equations is indeed crucial to understand the issue of the cosmological constant [1^. 

Finally we point out that the exponential potential in E frame always emerges even if 
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A in J frame is multiplied by a monomial 0^ (£ 7^ 4), whereas a power-law potential in 
E frame, as discussed in many references, can be interpreted as coming from a somewhat 
awkward choice of the J frame potential 0^(ln0)^ added to the prototype BD model. 
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Figure 1: The potential ^((T, $) given by (|^). Along the central valley with $ = 0, the 
potential reduces to the simpler behavior Ae~^''°" as given by (|l|), but with $ 7^ 0, it shows 
an oscillation in the a direction. 
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Figure 2: An example of the solution. Upper diagram: 6 = In a* (dotted), cr (solid) and 2$ 
(dashed) are plotted against logt*. The present epoch corresponds to logt* = 60.1 — 60.2, 
while the primordial nucleosynthesis must have taken place at logt* ~ 45. The parameters 
are A = 1,C = 1.5823, m = 4.75,7 = 0.8, k = 10, Cd = 0.005 (See Sect. 5 for the effect 
of the dust component). The initial values at ti — 10^^ are ai — 6.7544, a[ — {a, prime 
implies differentiation with respect to r = lnt^,),$i = 0.21, = — 0.002, pirad = 3.7352 x 
10^^^,Pidust = 4.0 X 10~^^. The dashed-dotted straight line represents the asymptote 
of a given by t/{2(). Notice long plateaus of a and $, and their rapid changes during 
relatively "short" periods. Middle diagram: a* = 6' = t^,H^, for an effective exponent in 
the local power-law expansion a* ~ t^* of the universe. Notable leveling-offs can be seen 
at 0.333, 0.5 and 0.667 corresponding to the epochs dominated by the kinetic terms of 
the scalar fields, the radiation matter and the dust matter, respectively. Lower diagram: 
log ps (solid), the total energy density of the cr — $ system, and logp*m (dashed), the 
matter energy density. Notice an interlacing pattern of ps and p*m, still obeying ~ 
as an overall behavior. Nearly flat plateaus of ps precede before ps overtakes p^rn, hence 
with ft\ passing through 0.5. 
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Figure 3: Magnified view of logp^ (solid) and logp*m (daslied) in tlie lower diagram of 
Fig. H around the present epoch. Note that ps is very flat in this diagram extending back 
to the past oi z = 5.2 — 6.9 for the assumed age (1.1 — 1.4) x lO^^^y. 
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Figure 4: Magnified view of a (solid) and 0.02$ + 44.25 (dashed) in the upper diagram 
of Fig. 0. Note that the vertical scale has been expanded by approximately 330 times as 
large compared with Fig. ^ With the time variable r = Int*, the potential V grows as the 
multiplying factor t1 = e^"^. The potential wall for a becomes increasingly steeper, thus 
confining a further to the bottom of the potential, as noticed by an oscillatory behavior 
of a with ever increasing frequency measured in r. (In this particular example, a was 
trapped finally long after it had been decelerated to settle near one of the maxima of the 
potential.) The growing V causes $ eventually to fall downward, resulting in the collapse 
of the confining potential wall. The stored energy is then released to unleash a. 
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Figure 5: An example of the solution without mini-inflation around the present epoch, 
to be compared with Fig. The qualitative change of the behavior is a consequence of 
choosing a slightly different value ai = 6.7764 as compared with 6.7544 in Fig. ^, with the 
rest of parameters remaining the same. In the upper diagram we flnd a making a smooth 
transition to the asymptotic behavior without being catapulted. Also $ approaches the 
value zero so slowly that it does not look it does in the limited range of time. 
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Figure 6: Magnified view of logp^ (solid) and Iogp*m (dasfied) in tlie lower diagram of 
Fig. ^ around the present epoch, to be compared with Fig. |^. No crossing is seen, with 
reaching only 0.29. 



15 



45.5 I r 




45.2 I ' ' ' ^ ' ' ' ' ' ' ' 1 

56 58 60 62 

log t 



Figure 7: Magnified view of a (solid) and 0.02$ + 45.25 (dashed) in tfie upper diagram 
of Fig. ^. Note that no oscillatory behavior is seen for a plotted in the same rate of 
magnification as in Fig. ^. This is due to the absence of the force which has protected a 
from flowing out as the potential grows as tl. 
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